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An iterface structure between turbulence and laminar flow is investigated in two-dimensional 
channel flow. This spatially localized structure not only sustains itself, but also converts laminar 
state into turbulence actively. In other words, this coherent structure has a functionality to gen¬ 
erate inhomogeneity by its inner dynamics. The dynamics of this functional coherent structure is 
isolated using the hltered simulation, and a physical perspective of its dynamics is summarized in a 
phenomenological model called an “ejection-jet” cycle, which includes multiscale interaction process. 
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INTRODUCTION 

Turbulence ubiquitously appears in nature: from 
qaurk-gluon plasma [T] to the Universe [2]. Because 
of its strong nonlinearity, most studies related to tur¬ 
bulence may have adopted more or less statistical or 
coarse graining approaches [3] . Though they have vividly 
revealed phenomenological and/or kinematic natures of 
turbulence such as the energy transfer among different 
scales and places, these statistical treatments are not suf¬ 
ficiently adequate to elucidate concrete mechanisms of 
even such fundamental processes of turbulence: For ex¬ 
ample, what substance, e.g. vortices, transfers energy or 
why the energy transfer occurs. On the other hand, the 
dynamical systems approaches to turbulence have helped 
us describe these mechanisms with numerically obtained 
components (invariant sets) in the phase space such as 
fixed points, periodic orbits and their connections [4]. 

Recent developments in the dynamical systems ap¬ 
proach to turbulence arrive at the next stage, where the 
spatial inhomogeneity is taken into account. Famous ac¬ 
tors on the previous stage are the “minimal” flows [5], 
which mean direct numerical simulations with minimal 
system sizes reproducing elementary processes and some 
statistical quantities of turbulence. The phase spaces em- 
beding them are effectively low-dimensional. However, 
those of the spatially inhomogeneous turbulent flows are 
no longer low-dimensional, and it is quite hard to treat 
such high-dimensional phase spaces both theoretically 
and numerically. 

One simple strategy to overcome this high- 
dimensionality is to consider spatially localized 
self-sustaining structures as building blocks (BBs) 
of turbulence. Though each of BBs may consist of 
internal fundamental elements, each block is expected 
to be effectively low-dimensional like the minimal flows. 
Indeed, various types of numerical exact solutions to 
the Navier-Stokes equation corresponding to localized 
coherent structures have been obtained so far in pipe 
flow [6] , plane Couette flow HH] and asymptotic suction 
boundary layer m- Then the dynamics of spatially 
extended systems are expected to be decomposed into 


that of each localized coherent structure and their 
interactions. 

At first glance this building block strategy may be in¬ 
compatible with global inhomogeneity since they are in¬ 
troduced to represent local dynamics. One scenario to 
treat global inhomogeneity in this strategy is to regard it 
as a collective dynamics among BBs. Since each of BBs is 
represented by a low-dimensional model, a coarse-grained 
model governing their interactions can be deduced as 
done for the chemical oscillations m- This phase re¬ 
duction m scenario has succeeded in explaining proper¬ 
ties of “puff” in pipe flow m- We demonstrate, however, 
that this scenario breaks down at least for a fundamental 
inhomogeneous system, namely turbulent-laminar inter¬ 
face. 

Instead, we try to deal with global inhomogeneity by 
extending roles of each BB. We focus on a turbulent- 
laminar interface in two-dimensional channel flow. As 
clarified in this paper, a localized self-sustaining struc¬ 
ture is embedded in the interface. This structure, which 
we call chaotic interface (Cl), produces turbulence down¬ 
stream by its inner dynamics while invading upstream 
laminar flow. The turbulent-laminar interface is gov¬ 
erned by Cl in this sense; i.e. this global inhomogene¬ 
ity is generated by the local dynamics. We introduce the 
term “functional” coherent structures (FCS) to represent 
such active localized coherent structures. We clarify the 
dyanamics of Cl in detail and evaluate how this func¬ 
tional building block scenario explains this global inho¬ 
mogeneity. 


CHAOTIC INTERFACE STRUCTURE 

Two-dimensional laminar channel flow has the same 
critical Reynolds number Rcc as three-dimensional one. 
In 2D case, the TS-wave solution appearing at this crit¬ 
ical point bifurcates into a weak chaotic state, which we 
call chaotic TS-wave, as its Reynolds number increases 
nMi]. In this paper, we consider a channel which con¬ 
tains turbulent downstream region and laminar upstream 
region. 
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FIG. 1. A snapshot of turbulent vorticity field. C varies from —25 to 25 on the walls, and regions C > 1-5 (C < —1-5) are 
colored by the same color of C = 1-5 (C = —1.5). 


We adopt a frame of reference moving at a speed c/ 
against the laboratory frame for Cl not to march. The 
streamwise and the wall-normal coordinates are denoted 
by X and respectively in this interface frame. The 
system is non-dimensionalized by the half width of the 
channel, so ^ G [—1,1]. u denotes the velocity field in this 
frame. We deal with a very long box [0, 207r] X [-1,1] 
periodic in x, to emulate the dynamics realized in an 
infinitely long channel. Since the walls move in the in¬ 
terface frame, the non-slip boundary conditions become 
±1) = —c/x, where x denotes the x directional unit 
vector. The Raynolds number Re is fixed to 8000 in this 
paper. 

To analyze the dynamics of this process in a finite com¬ 
putational box, we have to keep supplying laminar region 
since the turbulent region becomes wider as time goes on. 
We resolve this problem using the damping filter in 
the interface frame. We introduce a linear damping term 
into the incompressible Navier-Stokes (NS) equation to 
reproduce a laminar Poiseuille flow Ul = (1 — — ci)x 

in a small region Q = [0,1.4] X [-1,1]: 

d'lJL 1 

-^ + {u-V)u = -Vp + ^ {u - Ul) , 

where the last term of NS equation is the damping filter 
term. Since c/ = 0.855 is faster than the phase velocity 
of the chaotic TS-wave, this damping term laminarizes 
it, and the laminarized flow returns upstream due to the 
periodic boundary condition. 

In this setting a turbulent-laminar interface is simu¬ 
lated permanently. A snapshot is displayed in Fig. 
using the turbulent vorticity ( = (V x {u — This 

figure shows that there are three regions: weak turbu¬ 
lence {x ^ 20), chaotic interface (20 ^ x ^ 34), and 
laminar (34 ^ x) regions. Moreover, the chaotic inter¬ 
face contains dynamic inner structures. The snapshot 
shows a meandering bulk structure and strong wall shear 
layers. The weak turbulence consists of spatially mod¬ 
ulated chaotic TS-waves. The chaotic interface is noth¬ 
ing but FCS, and we will reveal in the following that 
it generates the weak turbulence. We first investigate 
the energy balance of these regions, and then construct 
a phenomenology for its self-sustaining mechanism and 
functionality. 

To focus on its streamwise inhomogeneity, we consider 
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FIG. 2. x-t plot of Pp, , and D^. The characteristic 
structure around x ~ 25 corresponds to the vortex ejection 
processes. 


the ^-averaged energy balance equation: 

OE 

+ dx{Ju + Jy) = Pp P Py — Dy + F‘ 

It should be noted that the energy is defined in the in¬ 
terface frame: E[xR) = j^-^dy\\u\\^/2. Since the walls 
move, there is an energy injection due to the viscosity on 
the walls Py = P+ pPy ^ where P^ = TcidyUx\y=±i/Re 
in addition to the bulk viscous dissipation 

-Di/ = ^ y dy(2 + {d^Uy)^ + {dyU^)'^^ . 

The term Pp{xR) = — f^^dy {u • V) p represents the en¬ 
ergy injection due to the pressure gradient, and takes 
both positive and negative values. Pp > 0 means the 
flow accelerated by the pressure gradient, and Pp < 0 
does the flow against the pressure gradient. Pp balances 
almost with the gradient of the energy flux dxJu and 
their spatial means are smaller than those of the viscous 
terms Py and Dy. The flux due to the viscosity Jy is 
negligible, and thus neglected hereafter. E is the energy 
damping by the filter term. The three terms Pp, P+, and 
Dy are displayed in Fig.|^ which illustrates the traveling 
of each structures. Reflecting the chaotic nature of the 
interface, these values are not exactly periodic. P^ (x, t) 
nearly equals to — P+(x, tpTpl2)^ where Tp ^ 15 denotes 
an approximate period of the recurrent motion at each 
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FIG. 3. A schematic view of the eject ion-jet cycle. 


point, and thus Pj^ is recurrent with the half period Tp/2 
like as Pp and Dj^. 

To confirm that the chaotic interface maintains itself 
in terms of energy balance, the energy balance equation 
is averaged over the interface region x G [20,34]: 

P 

1“ ^Ju ~ Pp,I T ^1^,1 

where the inferior •/ denotes the average over the in¬ 
terface, and AJu{t) = J^i(34,t) — Ju{20^t). The time 
average (•) of these terms are calculated: {dEijdt) = 
-6.7x10-^ ^0, {AJu) =3.0x10-^, (P^,/) =4.3x10“^, 
= 1-2 X 10“^, and {D^j) = 1.3 x 10“^. It should be 
noted that there is the averaged energy leak {AJy) > 0, 
which means that the chaotic interface is self-sustainable 
in terms of the time-averaged energy balance, and even 
an energy supplier to the weak turbulence. This energy 
leak reflects the functionality of the chaotic interface, i.e. 
the chaotic interface sustains the weak turbulence. The 
right after region of the chaotic interface has larger en¬ 
ergy or stronger turbulent intensity than the downstream 
side of the weak turbulent region or asympototic chaotic 
TS-wave. This convective relaxation process from this 
energy excess state to the asympototic chaotic TS-wave 
state is similar to a temporal relaxation process of a min¬ 
imal 2D channel flow, which is not shown in this paper. 
This similarity and the relationship to the phase reduc¬ 
tion scenario is left to future works. 


EJECTION-JET CYCLE 

Here we give a concrete description of the self- 
sustaining mechanism of the chaotic interface. This sus¬ 
taining process is constituted by the interaction among 
vortex ejections on the walls and the meandering jet in 
the bulk region. This collective dynamics is further split 
into three steps as summarized in Fig. In the step (i), 
a pair of sheet-like vortices is excited by the instability 
of the laminar flow near the wall triggered by the me¬ 
andering jet. The amplitude of the meandering decays, 
and the jet gets straight as going upstream. This sug¬ 
gests that a straight jet is convectively stable. Since the 
straight jet does not excite the vortex pair, it does not 
appear in x > 35. 

The step (ii) is the convective growth of the vortex pair. 
The thin vortex pair generated in the step (i) grows up 
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FIG. 4. The growth of the vortex pair. Three figures display 
different parts of the same snapshot. Regions C>3(C<—3) 
are colored by the same color of C = 3 (C = — 3) to emphasize 
the bulk vorticity, though the vorticity field y) reaches its 
maximum Cmax — 25 around x ~ 24. 


into an intense vortex ejection. This process is displayed 
in Fig. 1^ which picks up three continuing parts from a 
snapshot. Since the vortex pairs grow convectively, one 
snapshot of the entire channel gives three snapshots of the 
growing vortex pairs. In the energy viewpoint, it should 
be emphasized that this instability is not absolute but 
convective in both the interface frame and the labora¬ 
tory frame. In the laboratory frame, the vortex pair goes 
upstream at a constant speed Cy 0.5 > 0, and, in the in¬ 
terface frame, it goes downstream at Cy—Cj —0.35 < 0. 

The step (iii) is the vortex ejection process, which ex¬ 
cites the jet and makes it meander. The ejection oc¬ 
curs on the downstream side of the chaotic interface, 
namely around 22 < x < 28. Then the cycle is closed, 
and we call this cycle an “ejection-jet” cycle (EJC). A 
very strong shear accompanies this vortex ejection pro¬ 
cess. The wall unit ly is estimated at 2.1 x 10“^, and 
the friction Reynolds number Rcy = is about 460. 
This means that the width of the interface is 5000 times 
larger than ly. Therefore, we should regard this interface 
structure as a large-scale motion in the wall-turbulence 
context. After the intensive ejection process, the vortex 
structures are swept downstream, and this corresponds 
to the leak of the energy {AJy) from the interface to the 
weak turbulence region. 

To complete the EJC model, let us consider how the in¬ 
vading speed Cl is determined. There are two dynamical 
processes, the convective growth of the vortex pair and 
the decay of the jet meandering. Eirst, we suppose that 
both the traveling speed of each vortex pair Cy and the 
period Ty necessary to grow up are constant. Erom Eig.[^ 
we estimate them at Cy 0.5 ±0.05 and 20±2 [18]. 
Their inaccuracies are due to the inaccurate definitions of 
them, and more accurate and quantitative arguments are 
left to future works. Then d{ci) := \cy — ci\Ty denotes the 
distance between the birth point of the vortex pair and its 
ejection point. Next, we introduce a characteristic length 
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A of the decay of the meandering. Since this process is a 
nonlinear energy redistribution, we cannot define it from 
the spatial linear decay rate, but instead we measure the 
distance between the point where (max^ \uy{x^y)\) takes 
its maximum {x = 25.4) and the point where it becomes 
almost zero first {x ^ 38 ± 2). Since no vortex pair is 
excited when the jet does not meander, the EJC model 
requires these two length are equal: 

\Cy Cl\Ty = A. 

This condition connects two values characterizing the dif¬ 
ferent dynamical processes, and thus we should regard 
this condition as a self-consistent equation for c/. The 
above estimates are consistent with c/ = 0.855. 

Let us review the EJC model by introducing filtered 
simulations. We make other three runs in which the fil¬ 
tered region ft is set to damp one of the specific processes, 
namely, (a) weak turbulence, (b) vortex ejection, and (c) 
vortex pair excitation. Although we cannot split out each 
primary dynamics completely, these filtered simulations 
help us confirm the EJC model. These simulations use 
a snapshot of the previous simulation as an initial value, 
and animations visualized by the turbulent vorticity are 
included in the supplementary materials. 

Case (a): we set = [0,22] x [—1,1] to damp 

the weak turbulent region, and to confirm the self¬ 
sustainability of the chaotic interface. In this setting 
we yield a permanent chaotic interface, whose invading 
speed and the spatial structure are hardly changed. We 
conclude that the following weak turbulence is additional 
as assumed in the EJC model. Eurthermore, the selec¬ 
tion process of cj and the spatial structure is completely 
closed in the chaotic interface. In other words, the weak 
turbulence region plays no role in the selection process. 

Case (b): we set ft^^^ = [0,30] x [—1,1] to confirm 
that the jet is maintained by the acceleration due to the 
vortex ejection. If the meandering jet is self-sustaining, 
this simulation could yield a permanent finite amplitude 
solution. However, the laminar flow has occupied whole 
region. In this sense, the meandering of the jet is only 
a component mechanism of this self-sustaining process, 
and is not self-sustaining. 

Case (c): we set = [30, 207r] x [—1,1] to obstruct 
the step (i). In this case the non-filtered region of the 
chaotic interface (20 < x < 30) keeps alive on the same 
position until t < 20, and then it travels downstream. 
This time lag corresponds to the growth time Ty of the 
vortex ejection, and thus this result also supports the 
EJC model. After a long transient, another chaotic in¬ 
terface is reconstructed around 15 ^ < 27, and their 

invading speed and spatial structure are same as the pre¬ 
vious one. This result insists that the chaotic interface 
structure is robust while there is a laminar flow on its 
upstream. This robustness is an important issue for the 
pattern selection problem, but the current framework of 


the dynamical systems approach lacks tools applicable 
for settling the issue. 

CONCLUDING REMARKS 

We have investigated the self-sustainability and func¬ 
tionality of Cl in two-dimensional channel flow as an ex¬ 
ample of ECS, which yields the inhomogeneity between 
two asymptotic homogeneous states, upstream laminar 
flow and downstream chaotic TS-wave. We have intro¬ 
duced a phenomenology summarized in the EJC model, 
which consists of the vortex ejection and the meander¬ 
ing jet. The localized dynamics of Cl is isolated by the 
filtered simulation, and deconstructed by the energy bal¬ 
ance analysis. The EJC model well represents both the 
invading process on its front and its functionality in sus¬ 
taining the weak turbulence on its tail. As a result, how¬ 
ever, this functionality prevents us from obtaining an ex¬ 
act localized solution corresponding to the chaotic inter¬ 
face as done for various coherent structures because weak 
turbulence must attach to the interface. The damping fil¬ 
ter works effectively in isolating the localized dynamics 
of CL 

The self-sustaining mechanism described by the EJC 
model is also an example for collective dynamics of multi¬ 
scale structures. Different from Waleffe’s self-sustaining 
process m which utilizes an absolute instability, the 
EJC model does a convective instability, which needs 
a sufficient space to grow up. The convective instabil¬ 
ity makes it possible for the structures of different scales 
to interact with each other, namely the meandering jet 
of large scale and the wall shear of small scale. This 
multiscale interaction mechanism may be applied for the 
large-scale motion in three-dimensional wall-turbulence 
[20] , although the chaotic nature of Cl is far weaker than 
that of three-dimensional wall-turbulence. Eurthermore, 
it may also be a prototype for more general multiscale 
collective dynamics. 

We have introduced the functional coherent structure 
(ECS), which extends the well-known coherent structure 
perspective. Previous studies have focused on the self¬ 
sustainability of coherent structures, but we do on its 
additional functionality. We expect that the idea to as¬ 
sign functionalities of turbulence to localized coherent 
structures may work well for other cases. Energy and 
momentum transfers in fully-developed wall-turbulence 
are possible applications since functional Waleffe’s SSP, 
if it exists, may be embedded near the wall. Eor further 
development of the building block strategy, we will have 
to combine this functional building block scenario with 
the phase reduction scenario. In other words, we have to 
establish a framework involving phenomenological low¬ 
dimensional models of ECS and their interactions, and 
it is left to future works. This framework will be an es¬ 
sential tool for the dynamical systems approach to inho- 
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mogeneous turbulence and more general spatiotemporal 
chaotic systems. 
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